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Constraints of a pulsation frequency on stellar parameters in the eclips- 
ing spectroscopic binary system: V577 Oph 
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We present a preliminary spectroscopic analysis of the binary system V5770ph, observed during the summer of 2007 
on the 2.6m NOT telescope on La Palma. We have obtained time series spectroscopic observations, which show clear 
binary motion as well as radial velocity variations due to pulsation in the primary star. By modelling the radial velocities 
we determine a full orbital solution of the system, which yields Ma sin 3 i = 1.562 ± 0.012 Mg and Mb sin 3 i = 
1.461 ± 0.020 Mg. An estimate of inclination from photometry yields a primary mass of ~ 1.6 Mq. Using this derived 
mass, and the known pulsation frequency we can impose a lower limit of 1 Gyr on the age of the system, and constrain 
the parameters of the oscillation mode. We show that with further analysis of the spectra (extracting the atmospheric 
parameters), tighter constraints could be imposed on the age, metallicity and the mode parameters. This work emphasizes 
the power that a single pulsation frequency can have for constraining stellar parameters in an eclipsing binary system. 
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1 Introduction 

Large amplitude 5 Scuti stars are 1.5 - 2.5 M© mostly MS 
stars, with few observed pulsation modes. Interpreting these 
pulsations is difficult, mainly because 1) generally, the fun- 
damental parameters of the star are poorly constrained, and 
2) we do not know enough information about the mode char- 
acteristics (degree I, azimuthal m, and radial order n) to al- 
low us to use the mode as an extra observable constraint. To 
overcome these obstacles we can 1) study a pulsating star 
in a detached binary system, where the fundamental param- 
eters can be determined independent of using the pulsation 
mode, and 2) use new techniques to study the observational 
profiles of the oscillations to constrain its mode properties. 

V577 Oph is a detached eclipsing binary system with a S 
Scuti component, and is an ideal candidate for a seismolog- 
ical analysis. Its orbit is eccentric (e = 0.22), the dominant 
pulsation period of the 6 Scuti object is of the order of 1.3 
hours, l/100th the size of its orbital period, and the few pub- 



lished light curves indicate no out-of-eclipse modulations. 
All of this indicates that the system is detached, and so we 
may assume that the 8 Scuti evolves as a single star i.e. no 
mass transfer, tidal distortions or oscillations produced by a 
tidally-locked orbit. 

The system was first studied by Vol kov 19901 The au- 
thor concluded that the primary component (the more mas- 
sive and more luminous) is the pulsating star, due to the 
amplitude changes during primary and secondary eclipse. 
They also measured B-V = 0.488 mag outside of eclipse, 
0.504 mag during primary eclipse minimum, and 0.509 mag 
during secondary eclipse minimum. Using the photometric 
eclipse depths, the orbital period IT has been measured as 
6.079 days, and a dominant pulsation frequency of 0.0695 
days has been detected. The spectral types of both compo- 
nents have been estimated as A8 from interstellar reden- 
ning (Shugarov 1985 — referenced in Diethelm 1993), and 
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more recently, Zhou (2001) obtained new photometric light 
curves and extracted the dominant pulsation frequency of 
14.390 cycles per day with an oscillation amplitude of 0.029 
mag. 
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In this paper, we discuss our recent spectroscopic obser- 
vations of V5770ph. As an initial analysis, we extract the 
the radial velocities and determine the orbital parameters. 
Using an estimate of inclination from photometry to obtain 
the mass, we use seismic models to discuss the global prop- 
erties of V577 0ph. 



2 Observations 

We obtained 5 consecutive nights on the Nordic Optical 
Telescope (NOT), located in the Observatorio del Roque 
de los Muchachos on La Palma, Canary Islands. The NOT 
is a 2.6m telescope with an alt-azimuth mount. The instru- 
ment used is the Fibre-fed Echelle Spectrograph (FIES) — a 
cross-dispersed high-resolution echelle spectrograph with a 
maximum resolution of R=65,000. Our setup was medium- 
resolution with R=45,000, with spectral range: 3640 - 7455. 
The maximum system efficiency (detector-instrument-fiber- 
telescope-atmosphere at airmass 1.0) is 9% at 6000 A. 

During the nights of July 4 - July 8, 2007 inclusive, 
we obtained spectra of V577 Oph, using 10 minute expo- 
sures. Accounting for calibration, Targets of Opportunities 
and other interruptions, each night an average of 32 science 
images in sets of 4 for V577 Oph were taken. Due to the 
faintness of our source ( V = 1 1 .0) and the potential of over- 
exposing ThAr, we opted to obtain the calibration frames 
at regular intervals of about 40-45 minutes. These frames 
were used to calibrate the wavelength of the two images be- 
fore and after this frame. This results in a loss of accuracy in 
wavelength, but of the order of 0.002 A at A = 3659.6294 A. 
This corresponds approximately to an error of 0. 1 - 0.2 km 
s- 1 . 

We used the FIEStool (Stempels, 2004)3 reduction 
software to calibrate the wavelength and reduce the spectra. 
This tool is optimized for FIES data, but the reduction is 
standard and calls I RAF to perform some of the tasks. We 
subsequently shifted the spectra to barycentric frame. 



3 Analysis 

3.1 Determining the radial velocities 

Because the spectral lines are not very broad, and the bi- 
nary motion can be easily seen, we chose to measure the 
radial velocities by modelling a set of absorption profiles 
with Gaussian functions, where the central position of the 
Gaussian corresponds to the measured radial velocity. We 
identify the lines by eye, by inspecting the time series such 
as that shown in Fig. Q] We chose the lines that are not 
blended with other spectral features e.g. those around 5535 
A in Fig. Q] The lines between 5525 and 5530 A were not 
chosen, because near relative radial velocities of zero, it is 
difficult to measure these blended lines. 



o 




JZ 
*Q0 




5520 



5535 
-10„ 



5540 



5525 5530 

Wavelength (10~ lu m) 

Fig. 1 Small spectral range showing the evolution of the 
doubled spectral lines. In some places it is clear to see the 
evolution of one line, while in others the lines are blended 
especially at maximum separation. 

Figure |2] (top panel) shows the fitted radial velocities 
using the spectra in a range around 5183 A as a function of 
phase cj) from an arbitrary reference phase. The fitted values 
are shown within the dashed lines, which represent the be- 
ginning and end of one orbital period. At <f> < and <fi > 1, 
the same data are plotted modulo II (orbital period). The 
large noise associated with the primary component (posi- 
tive radial velocities for <fr < 0.4) is the dominant oscilla- 
tion mode with an amplitude of about 5 km s _1 . An initial 
inspection of these results indicates that the components of 
the system are of similar mass — the spectral features also 
show evidence of this. 

3.2 Orbital solution 

We determined the orbital solution by fitting simultaneously 
the free parameters: the semi-amplitude of the radial veloc- 
ity curves for component A and B, Ka and Kb, the sys- 
temic velocity 7, the eccentricity e with argument of peri- 
astron u>, as a function of either time ti or phase <f>i, where 
<f>i = t ' m ° dn - T and T — is an arbitrary reference 
value. When we fitted as a function of phase, the reference 
phase 0o was fitted with II fixed, and when fitting with time 
we fitted both T and IT. 

We fitted the data using a \ 2 te st: 
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Table 1 Orbital Elements from Radial Velocity Curves 







Po 


Pab 


Pa 


Pb 




Pab 


Pa 


Pb 
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e 




0.2 


0.212 


0.169 


0.223 


/ f\ f\f\ A \ 

(0.004) 


0.204 


0.167 


0.224 


(0.006) 


0.204 


7 


(kms" 1 ) 


-37.0 


-38.38 


-38.26 


-37.62 


(1.34) 


-38.24 


-38.04 


-37.66 


(0.87) 


-38.24 


CO 


(°) 


47.0 


50.13 


46.39 


47.37 


(0.23) 


47.80 


47.59 


47.69 


(0.22) 


47.80 


K A 


(kms^ 1 ) 


85.0 


83.78 


83.53 


83.79 


(0.60) 


83.99 


83.49 


83.99 


(0.58) 


83.99 


K b 


(km s _1 ) 


90.0 


89.65 


89.65 


90.05 


(0.27) 


89.79 


89.79 


90.06 


(0.35) 


89.79 


To 


(days) 


2.0 


1.727 


1.669 


1.702 


(0.017) 


1.699 


1.682 


1.706 


(0.01) 


1.699 


n 


(days) 


6.079 


6.145 


6.053 


6.072 


(0.030) 










6.079 


Xr 






0.649 


0.691 


0.523 




0.672 


0.687 


0.520 







where Oj, ej and Mi are respectively the N observations, 
measurement errors and the model observables produced 
from a specific set of model parameters. An initial guess of 
the parameters are needed, and these are obtained by using 
the published information about the photometrically derived 
II and e, and directly inspecting the radial velocity curves. 
The initial values are shown in the first column of Table Q] 
The convergence criteria was given by setting Xr (reduced 
X 2 ) to 1 as well as testing how the free parameters changed 
with each iteration of the minimization algorithm. When the 
parameters were considered "stable" the minimization was 
usually stopped. 

The fitting was performed by modelling the individual 
radial velocity curves i.e. just one of Ka and Kb and by 
modelling both simultaneously. These are denoted respec- 
tively in Table [T] as the column headings Pa, Pb and Pab- 
Note that when modelling only one curve, and when the pa- 
rameters converged, we fixed this solution to subsequently 
determine Ka or Kb from the data of the other compo- 
nent. The first set of columns show the results when we fit- 
ted II and To as a function of time, while the second set of 
columns show the results when we use phase (note that the 
parameter To has the meaning of <f>o for this second set of 
solutions). The fourth column lists the average value of the 
uncertainties from these fits. 

4 Results 

The fitted parameters and uncertainties (<7j) are shown in 
Table Q] (see above paragraph for details). 7, u>, Ka and 
Kb return stable values while using all six fitting functions 
(data columns 2-4 and 6 - 8): all values fall to within 
their quoted uncertainties. Additionally there is no qualita- 
tive difference between these fitted parameters while leav- 
ing II free or fixing it, the only variation found is the value 
of the uncertainty on 7. 

The initial value of eccentricity was estimated as e = 
0.22 (Shugarov 1985). When we fitted the radial velocities 
using the individual curves and then using both together we 
obtain very different values: 0.207, 0.168 and 0.224, with 
these values almost reproduced while fitting with II fixed. 
Figure [2] shows the fitted radial velocity curves for each of 
the second three solutions. The green, red and blue curves 



Table 2 Physical Stellar & Orbital Parameters 





P 


a 


M A siir* i (M©) 


1.562 


0.012 


M B sin 3 i (M Q ) 


1.461 


0.020 


asini (AU) 


0.0942 


0.0003 


II (days) 


6.0791 


0.0001 



show the orbital solutions from the fits using both sets of ra- 
dial velocity data, the primary component only and the sec- 
ondary component only, respectively. To account for the ob- 
served dominant oscillation amplitude, 5 kms -1 was added 
to the measurement errors for the primary component for 
the fitting process. 

Obtaining different fitted values for IT will affect the fit- 
ted value of To. This is the reason that a larger difference 
between the To values is seen for the first three solutions, 
than for the second three solutions. While keeping II fixed, 
the variation in To is only within 2a. 

Because of the large amplitude residual, it would be less 
feasible to regard the fit using the pulsating component data 
only as the optimal (red curve). Inspecting the data points 
(for the secondary component) and the red curve indicates 
that the fit is not optimal for both components. For this rea- 
son, the choice of e = 0.16 is discarded. 

The fitted values while using the data from both stars 
is adopted (green curve), and are given in Table Q] last col- 
umn under "P". This solution is consistent with the solution 
while fitting just the data from the non-pulsating compo- 
nent (blue). Additionally, there is no difference in the de- 
rived values of a sin i, Ma sin 3 i and Mb sin 3 i when we fit 
II and when it is fixed. The lower panel of Fig. [2]shows the 
residual for the data from the primary pulsating component 
divided by the observational errors (including the additional 
5 km s _1 added to the observed error) in red, and for the 
secondary component in blue. All of the data points fall to 
within 2e indicating an adequate fit. 

Using the orbital solution given under the column head- 
ing P from Table [T] we determine the system's physical pa- 
rameters: a sin i = 0.0935 ± 0.0003 AU, M T sin 3 i = 3.024 
± 0.023 Mq, and q = M B /M A = 0.939 ± 0.006. Tabled 
summarizes the component and orbital parameters. 
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Fig. 2 Top: The orbital solution from fitting the radial 
velocity data. The data within the vertical dashed lines are 
the true data points, those at phase <f> < and <f> > 1 are 
the data repeated modulo II. The green, red and blue curves 
represent the resulting fit curves from Pb, Pa and Pab- 
Bottom: Residual of data minus fit radial velocity curve with 
Pab for the primary (red) and the secondary (blue). 



5 Seismic inference 



An estimate of the inclination of the system from photome- 
try (i > 80°) yields a primary mass of ~ 1.6 M (Zhou, 
priv. comm.). Figure [3] show theoretical frequencies of a 
1.63 M Q star for 2 values of metallicity: Z (initial metal 
mass fraction) = 0.020 and 0.025. These mode frequencies 
were calculated using the ASTEC and ADIPLS stellar struc- 
ture, evolution and pulsation codes (Christensen-Dalsgaard 
2008a,b). The lines correspond to fundamental (continuous) 
and first overtone (dashed) oscillation modes of different de- 
gree {I = 0,1,2 from bottom to top). The boxes represent 
the observed value of the frequency which includes an error 
to take into account that the frequency could have a non- 
zero azimuthal order (m where — I < m < I). 

If we assume that no g-modes are present, since these 
are not common among large amplitude S Scuti stars, Figure 
3 then shows that the system must have a lower age limit 
of 1 Gyr (younger stars exhibit higher frequencies). If we 
knew which mode was observed we could further constrain 
the age to either 1.2 - 1.3 Gyr or 1.6 - 2.0 Gyr. Also we 
see that there are three possible solutions to the identity of 
the mode, without using any observational technique to do 
this. These possibilities are: radial mode (fundamental or 
first overtone) or fundamental non-radial mode 1=1. We 
can also see that a metallicity measurement and a mode- 
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Fig. 3 Frequency values for various low degree modes for 
a 1 .63 M0 star (V577 Oph A) as a function of age for Z = 
0.020 (top) and 0.025 (bottom). See text for details. 

identification would provide further constraints on the age 
and chemical composition. 

6 Conclusions 

We obtained 5 nights of spectroscopic observations of the 
eclipsing binary with a S Scuti component V577 Oph. We 
extracted the radial velocities of both components, which 
shows very clear orbital motion of both stars and pulsa- 
tions in the primary (more massive) star. We fitted the ve- 
locities to solve for the orbital parameters and then derive 
for the components M A sin 3 i = 1.562 ± 0.012 M and 
M B sin 3 i = 1.461 ± 0.020 M©. An estimate of inclination 
implies a primary pulsating component of ~ 1.6Mq. 

A dominant oscillation frequency has been confirmed, 
and comparing this with some stellar evolution and pulsa- 
tion models allows us to impose a lower limit of 1 Gyr on 
the age of the star, and constrains the identification of the 
mode to either radial mode (fundamental or first overtone) 
or fundamental non-radial mode I = 1. We plan to disen- 
tangle the spectra in order to determine the atmospheric pa- 
rameters of T ff, [M/H], and log <?. This information would 
additionally constrain both the age of the star and the pulsa- 
tion mode parameters. 
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